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We study the interference between the Fano and Kondo effects in a side-coupled double-quantum- 
dot system where one of the quantum dots couples to conduction electron bath while the other 
dot only side-couples to the first dot via antiferromagnetic (AF) spin exchange coupling. We apply 
both the perturbative renormalization group (RG) and numerical renormalization group (NRG) 
approaches to study the effect of AF coupling on the Fano lineshape in the conduction leads. With 
particle-hole symmetry, the AF exchange coupling competes with the Kondo effect and leads to a 
local spin-singlet ground state for arbitrary small coupling, so called "two-stage Kondo effect". As 
a result, via NRG we find the spectral properties of the Fano lineshape in the tunneling density of 
states (TDOS) Pc(i^) of conduction electron leads shows double dip-peak features at the energy scale 
around the Kondo temperature and the one much below it, corresponding to the two-stage Kondo 
effect; it also shows an universal scaling behavior at very low energies. We find the qualitative 
agreement between the NRG and the perturbative RG approach. Relevance of our work to the 
experiments is discussed. 

PACS numbers: 72.15.Qm,7.23.-b,03.65.Yz 



I. Introduction 

Fano resonance is the quantum interference effect be- 
tween a localized state with finite-width and a conduc- 
tion band 1]. The hallmark of the Fano resonance is 
the asymmetric lineshape in tunnelling density of states 
(TDOS) of the conduction band. One example of Fano 
resonance is the transport through low dimensional elec- 
tronic (Fermi) system with local impurities. The Kondo 
effect [4] plays an important role if these impurities carry 
unpaired spins. Recently, there has been growing interest 
both theoretically and experimentally in the Fano reso- 
nance associated with the Kondo effect via the STM mea- 
surements of noble metal surfaces [^-Q as well as in quan- 
tum dot devices [lol 11 1. The Fano resonance in these 



systems in general arises from two quantum interference 
effects: 1. between the broadened local level and the 
continuum conduction band and 2. between the Kondo 
resonance in the local level and the conduction band. 
The combined two effects give rise to rather complicated 
lineshpe in STM measurement of the TDOS. The Fano 
resonance in TDOS of conduction electrons in such sys- 
tems can be served as an alternative approach to study 
the Kondo effect in addition to the local density of states 
of the quantum dot. The Fano lineshape in TDOS of 
conduction electrons in the leads of a single Kondo dot 
system has been extensively studied, and it is sensitive 
to both the spatial phase of the free conduction electrons 
and the scattering phase shift associated with the Kondo 
effect. 

Very recently, the Fano resonance has been extended 
experimentally and theoretically [l^, [3 the side- 
coupled double quantum dot system where the com- 



petition between Kondo and Fano effects gives rise to 
change in conductance profile. In this paper, we inves- 
tigate the Fano-Kondo interference in the side-coupled 
double-quantum-dot systems where only one of the two 
dots (dot 1) connects to the leads while the other iso- 
lated dot (dot 2) is side-coupled to dot 1 16|, ll7|. In the 
Kondo limit where charging energy on each dot is large, 
an antiferromagnetic (AF) spin exchange (RKKY) cou- 
pling is generated via the second-order hoping between 
two dots competes with the Kondo effect, leading to lo- 
cal spin-singlet ground state for arbitrary finite values of 
J, so-called "two-stage Kondo effect" 3) 121- Previous 
studies on the side-coupled double-dot systems have been 
mostly focused on the dip of LDOS on dot 1 upon ap- 
plying the AF RKKY coupling. However, little is known 
about the feedback effect of the two-stage Kondo effect 
mentioned above on the TDOS of conduction electrons 
in the leads. In this paper, we generalize the Fano line- 
shape in TDOS of electrons in the leads as a result of the 
two-stage Kondo effect in side-coupled double-quantum- 
dot system. The systematic perturbative and numerical 
renormalization group approaches are applied here in the 
cases both with and without particle-hole symmetry. We 
find as a consequence of the two-stage Kondo effect, the 
spectral property of the Fano lineshape in TDOS of the 
leads develops an asymmetrical double dip/peak struc- 
ture; it also shows an universal scaling behavior at very 
low energies. We compare our NRG results with the per- 
turbative RG analysis. 

II. The Model Hamiltonian. 

Our starting Hamiltonian for the side-coupled double-dot 
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system is the single-impurity Anderson model for dot 
1 with additional antiferromagnetic spin-exchange cou- 
pling between dot 1 and the isolated dot 2 which side- 
coupled to it[l6|. 
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where and Ir denote the tunneling amplitudes to the 
left and right leads, respectively, and cj^^^ creates an elec- 
tron in lead a = L, R with spin a. This tunnel coupling 
leads to a broadening of the level on dot 1, the width of 
which is given by F = + Tj^ = 2i:{t\Q]^ + t'j^QR), with 
Ql/r the density of states in the leads. Here, i = 1,2 

labels the two dots, and Si = {^/2)J2aa' 'Aa'^crcr'd^^, is 
their spin. Each dot is subject to a charging energy, 
Ui ^ U2 — U — Ec- In the presence of particle- hole 
symmetry, we have edi — Note that in the Kondo 

limit where the charging energy Ec is large, the direct 
hoping between the two dots are strongly suppressed and 
an antiferromagnetic spin exchange coupling J > is 
generated via the second-order hoping processes. 

The physical observables of our interest are: (i). LDOS 
of impurity on dot 1: pdii^jj) = -^Ii^Gdii'^) ^'^'^ (ii) 
the TDOS of the conduction electron pdi^)'- Pc{^) = 
Pq + Spc{(^) where po is the density of states of the 
bare conduction electron leads : po = — /mG2(a; — 0) 
with ^^{u! — irj) being the bare conduction electron 
Green's function, and Spc{io) is the correction to the 
LDOS of the conduction electron due to the coupling 
between leads and the quantum dot system: Spduj) ~ 
^ImSGc{u! — irj). Here, the correction to the conduc- 
tion electron Green's function 6Gc(oj — irf) is given by: 

SGAi^-ir^) = —G°(u}~ir])Gdi(oj-i7])G"(uj-i7]) (2) 

Using Eq. [21 we haveQ 

Spc{oj) = ~Tpo X 

[{ql - l)ImGdi{oJ - iTi) - 2q,ReGdi{uj - zt?)], (3) 



with Qc being defined as 



ReG'liuj - ir]) 
ImG^(uj — i-q) ' 



(4) 



and it can be treated approximately as an frequency- 
independent constant (3;, .4]. Following Ref. [ill, below 
we apply both perturbative renormalization group (RG) 
and numerical renormalization group (NRG) approaches 
to calculate these quantities in the presence of particle- 
hole symmetry. Though the LDOS on dot 1 (or equiva- 
lently the imaginary part of the Green's function on dot 



1, ImGdi{u!)) at finite RKKY coupling J via both RG 
and NRG has been computed in Ref. [16[, the real part 
of Gdiiuj), ReGdi{(jj), which is also needed to analyze the 
spectral property of the Fano lineshape in the TDOS of 
the conduction electron leads {pc{uj)), has not yet been 
calculated by either perturbative RG or NRG approach. 
In the following, we provide a numerical and analytical 
analysis on the Fano lineshape for pd^o) by analyzing 
both the real and the imaginary parts of Gdi (w) at finite 
J via NRG and compare them with those via perturba- 
tive RG approach. 

First, we discuss the case for J — 0. For J = and in 
the presence of particle-hole symmetry (edi — ~Ui/2), it 
has been known that in the Kondo regime uj <^ T^. with 
Tk ~ Dqc^^^^^^ being the Kondo temperature for dot 1, 
Gdi{ijj) is well approximated by the single Lorentzian[l6|: 



Gdi(w) 
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with z — c ^ being the quasi-particle weight at the 
Fermi energy, and Tk — zT — c Tk being an energy of 
the order of the Kondo temperature, Tk- The precise 
value of the universal constant c relating Tk and Tk de- 
pends on the definition of Tk- Here, we define Tk as 
the half- width of the transmission T{uj) = —TlmGdiioj)- 
From fitting Gdi{uj) with the NRG data, we get c « 0.5. 
Note that by Fermi-liquid theory and principles of renor- 
malization group, the RKKY interaction also gets renor- 
malized by the same z— factor: J — >■ 7 16|. Here, J is 
slightly different from J due to the large logarithmic tail 
in ImGdi{i^)- The value of J is obtained from the fit 
of ImGdiico) to NRG data: J « l.lj[i6|. However, for 
uj >Tk, the above simple Lorentzian approximation fails 
to account for the large logarithmic tail in ImGdi{oj)- 
Therefore, corrections to the single Lorentzian approx- 
imation are needed in this case to more accurately de- 
scribe Gdi{to). Via the Dyson equation approach, tak- 
ing into account the interference between the Kondo res- 
onance and the broadened impurity level, we obtain a 
more accurate description for the Green's function of the 
dot iH: 

Gdi {lo) = G°i {uj) + G°i {Lo)fdi (c.)G°i {lo) (6) 

where the bare Green's function on dot 1, G^]^(w), de- 
scribing a local impurity level with a level broadening F 
and LDOS pda = — 0), is given by: 

l-n/2 , n/2 



^dl — 



LO — Edi + iV LO — Edi —Ui + iV 



with n —< n^-^ + n^j^ > being the average occupation 
number on dot 1; and Tdi{Lo) is the scattering T— matrix 
corresponding to the Kondo resonance, given approxi- 
mately byQ: 



Tdi{uj - irf) 



be 



i2S 



EK + i Tk + irj 



(8) 
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with b being a fitting parameter to be fitted with the 
NRG data for ImGdi{uj). In the presence of particle- 
hole symmetry, we have n — I, — 0. Here, d in Eq. [5] 
corresponds to the phase shift associated with the Kondo 
resonance scattering, and it gives 5 = 7r/2 in the case of 
particle-hole symmetry. By fitting Eq. IH] with the NRG 
data, we find b w z/{Trpdo)'^, which is in good agreement 
with the known result: —ImGdi{ui = 0) = l/F for a sin- 
gle impurity Anderson model[2]. In the Kondo regime 
(w <^ Tk and Ec » F) of our system and for J = 0, 
the bare Green's function on dot 1, G^^{uj), are approx- 
imately given by: ReG°^j^{io) « 0, '-^""^^^^"^ « pdo- The 
above approximations lead to the following approximated 
expressions for Gdi{i^) after including the interference be- 
tween the Kondo resonance and the broadened impurity 
(dot 1) level via Eq. [SI 



ReGdi{uj) w {TTpdoY 



buj 



ImGdiiuj) w -Trp^jo - (vrprfo) 



2 bTK 



(9) 



with -i/TOG°i(w = 0) being the LDOS of dot 1 

at w = 0. From Eq. [3] and Eq. [51 in the Kondo limit the 
correction to conduction electron density of states can 
therefore be expressed in terms of the well-known Fano 
lineshapefi, 0]: 



5pc{^, J 0) « po(- 



2g,e - 1 



1 



+ /3), (10) 



where e 



Tk 



,(3 = npdo^iql — 1). Note that in gen- 



eral the Dyson equation approach in Eq. [6] is also valid 
for both u! ^ Tfc and lu Tk in the presence of large 
particle-hole asymmetry: \edi — epl < F (with ep being 
the Fermi energy of the leads) where the the interference 
between the Kondo and broadened impurity level plays 
an important role in Gdi(cj)[3|. 



III. Perturbative Renormalization Group 
analysis 

Now, we turn on a finite RKKY coupling J. Follow- 
ing Ref. [3], to gain an analytical understanding we 
employing the perturbative renormalization group anal- 
ysis in the limit of J — t- 0. We restrict ourselves the 
case with particle-hole symmetry. Though some of the 
aspects in this case has been studied in Ref. [3], it 
proves to be useful to summarize its key results for fur- 
ther calculations on the Fano lineshape for pd^jj) in the 
presence of RKKY coupling J. In the limit J — >■ 0, "two- 
stage Kondo screening" takes place[il[i3|: The spin of 
dot 1 first gets Kondo screened below Kondo tempera- 
ture Tk ~ D e~'^^/^ with D being the bandwidth cut- 
off, the first stage Kondo effect. Then for energy scale 




FIG. 1: (a). ImGdi(ijj) (normalized to ImGdii^jJ — 0,J = 
0) « — 1/r) of dot 1 with particle-hole symmetry for dif- 
ferent antiferromagnetic RKKY couplings J calculated by 
NRG (solid lines) and perturbative RG (dashed lines) via 
Eqs. [6l and [TSl (dashed lines). Dotted lines are RG fits via 
Eqs. [5l and [TSl The NRG parameters are Ui = U2 = Do, 
Edi = ed2 = -0.5Do, r = 0.2L>o with Do = 1. For 
J = 0, we find Tk « 0.005-Do. The fitting parameters 
c « 0.5, and J ~ 1.1. (b). ImGdi{i^/T*) (normalized to 
ImGdiis^ = 0, J = 0) ~ ^1/r) shows an universal scal- 
ing behavior for ui < T* . The dot-dash line is the power- 
law {uj/T'f fit to the crossover function oi ImGdiij-^ /T*) for 
cj < T*, see Eq. (HI 



much below T^, the second stage Kondo effect occurs at 
uj < T* ^ Tk between dot 1 and 2 via the antiferro- 
magnetic RKKY coupling J where the spin on the dot 
2 gets Kondo screened. Here, the Kondo resonance peak 
in electron density of states on dot 1 plays the effective 
fermionic bath for the second stage Kondo effect. We will 
discuss how the Fano lineshape for Pc{uj) is affected in the 
presence of the antiferromagnetic RKKY coupling. Sum- 
ming up all leading logarithmic vertex diagrams leads to 
the followin g sc aling equation for the dimensionless ver- 
tex function [16| 



d(7(^,Tfe)) _ digiL^)j) 



dl 



dl 



(11) 



with the scaling variable defined as I = \og{TK /T'j^). 
Here, g{Lo) = g{oj)/z — ^ImGdi{ijj) is the rescaled ef- 
fective density of states of dot 1. Integrating this dif- 
ferential equation up to / = log(TR-/w), one obtains the 
dimensionless vertex function in the leading logarithmic 
approximation [lij : 



7(w,Tk) 



1 



|ilog^+logM 



with the second scale T* defined as 



T* = Tk exp(-^ Tk/J) 



(12) 



(13) 
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FIG. 2: ReGdi{i^) (normalized to -ImGdi{oj = 0, J = 0)) of 
dot 1 for different antiferromagnetic RKKY couplings J by 
NRG (solid lines). The dashed line is a fit to the NRG data 
for J = Q via Eq. [6] The other parameters are the same as in 

Fig.m 
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J=0.001 RG fit via Eqs. (5), (15) 
J=0.0015 RG fit via Eqs. (5), (15) 
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FIG. 3: ReGdi{i^) (normalized to -ImGdi{i^ = 0,J = 0)) 
on a logarithmic scale of uj/Tk for different antiferromagnetic 
RKKY couplings J by NRG (solid lines) and perturbative RG 
via Eqs. l6l and [151 (dashed lines). Dotted lines are RG fits via 
Eqs. [51 and [T^ The other parameters are the same as in Fi^ 



The second order self-energy correction to the retarded 
Green's function G^j simply gives the expression [l6| 



(14) 



where S=l/2. The Green's function of dot 1 after includ- 
ing self-energy and vertex correction is given bv[l6j|: 



(w) "--^^ '-Gii{uj) 



(15) 



where J(w) is replaced by ^[uj)/ q{uj)^ and Gdi{^jj) is given 
by either Eq. [5] (the Dyson equation approach) or Eq. [6] 
(the single Lorentzian approximation). Note that due to 
the logarithmic corrections in 7, the spectral density of 
dot 1 develops a dip at energies w ~ T* <g; Tk for any 
infinitesmall J, which suppresses the low-energy trans- 
mission coefficient through dot 1. Physically, this comes 
from as a consequence of the fact that electrons of energy 
ui < T* are not energetic enough to break up the local 
spin singlet and therefore their transport is suppressed. 
For a finite RKKY coupling J > 0, the real and imagi- 
nary parts of G'^i{uj) obtained in Eq. [15] via perturbative 
RG approach lead to an analytical expression for the cor- 
rection to the LDOS on dot 1, Sp'^{u!): 

6pi{u) = -rpo[{ql-l)ImGi,{u)-2qMeGi,{w)]. (16) 

Below we present the results via NRG with fits by the 
perturbative RG calculations. 



IV. Comparison to the Numerical 
Renormalization Group (NRG) analysis. 




FIG. 4: ReGdiioj) (normalized to -ImGdiioj = 0,J = 0)) 
on a logarithmic scale of uj/T* for different antiferromagnetic 
RKKY couplings J by NRG. The other parameters are the 
same as in Fig. [1] The dot-dash lines are power-law {uj/T*) 
fits to the universal crossover function of ReGdi{oj) for ui <^ 
T*, see Eq. [19] 



We have performed the NRG calculations on the sys- 
tem in the presence of particle-hole symmetry. The 
NRG parameters we used are: Ui = U2 = Dq — 1, 
edi = ed2 - -0.5, A = 2, Fl = r^. = 0.1 with Da being 
the bandwidth of the conduction electron baths. (Here, 
we set Do = 1 as the unit of all parameters.) Within 
each NRG iteration, we keep the lowest 1000 states. For 
J = 0, we find Tk « O.OOSDq- As RKKY coupling is 
increased, both real and imaginary parts of Gdi{i^) get 



splited at Ld = 0. First, as shown in Ref. [16|, the 
imaginary part of Gdi{uj) (proportional to DOS of dot 
1) at finite J shows a dip below the characteristic en- 
ergy scale T* for any arbitrary J > (see FiglIJa)). For 
small RKKY coupling J, the NRG results for ImGdi{io) 
can be fitted reasonably well by the perturbative RG ap- 
proach over an intermediate energy range T* <^ uj <^Tk. 
Furthermore, a clear universal scaling behavior of the 
KT type is observed from the NRG results of ImGdi{oj) 
for a; < T*: ImGdi{cu) « gogi^j/T*) (see Fig. [^b)) 
|16| . With particle-hole symmetry, the scaling function 
q(uj /T*) is completely universal. As pointed out in Ref. 
|16| , the ground state of the system at any finite J is a 
local spin-singlet (a Fermi liquid), the very low energy 
crossover of ImGdi{<^) for uj <^Tk vanishes as (cj/T*)^, 
following the Fermi liquid behavior: 




TImGdi{cd) 



ai(^) 



(17) 



FIG. 5: The Fano lineshape for pc(i^) (in unit of po) for dif- 
ferent antiferromagnetic RKKY couplings J by NRG. The 
dashed line is a fit to the Fano lineshape form Eq. [10] for 
J=0. The other parameters are the same as in Fig. 1. 



where a-\ k, 3.0 from the fit to the NRG data (see Fig. 
H (b)) [16|. Note that we find the perturbative RG 
approach via Eq. [HI leads to a better fit to the NRG 
results for ImGdi{i^) than that via Eq. [5l as expected. 

We now discuss the real part of Gdi{'-^)- For J = 0, 
ReGdi{<^) is antisymmetric with respect to w = and 
it shows a peak/dip at w « ±Tfc, signature of the first 
Kondo effect. As the RKKY coupling J is increased, 
the magnitude of the peak/dip in ReGdi{<^) for a; w 
decreases, indicating the Kondo effect is suppressed. At 
a much lower energy scale, T* Ri uj Tk, the Kondo 
dip-peak structure in i?eGdi(w) gets a further split with 
a width D w 2T*: it develops a negative- valued dip for 
Lj K, T*] while it shows a positive- valued peak for w w 
—T*. In the w — ^ limit, both positive and negative 
branches of ReGdi{uj) vanish (see Fig. [2] and Fig. [3]). 
We can get an analytical understanding of this behavior 
as follows: In the Kondo regime uj <^Tk, the real part of 
G'^^{ijj) is approximately given by (see Eq. [T5|) 



ReGi,{u) 



zuj{l 



16T? ' 



W2(l 



16T? 



(18) 



Prom the perturbative RG results, as w T*, J{uj) 
diverges, leading to the vanishing LDOS. As lo decreases 
to T* from above, the factor 1 - in Eq. [H first 



becomes negative then it approaches as w further ap- 
proaches 0. This explains the additional dip-peak struc- 
ture seen for T* in the NRG results. This qual- 
itative feature can be captured by the perturbative RG 
approach. However, the magnitudes of the dip-peak fea- 
tures via perturbative RG are much smaller than those 



J=0.001 RGfitviaEqs. (5), (15) 
J=0.0015 RGfitviaEqs. (5), (15) 
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FIG. 6: The Fano Hneshape for pc{(jj) (in unit of po) on a log- 
arithmic scale of w/Tfe for different antiferromagnetic RKKY 
couplings J by NRG (solid lines) and perturbative RG via 
Eqs. [6l and [T5l (dashed Unes). The dot-dashed lines are fits 
to the Fano hneshape form via Eq. [in] for J = 0. The dotted 
lines are the RG fits via Eqs. [5l and [151 The other parameters 
are the same as in Fig. 1. 



obtained from NRG. We believe the reasons for the devi- 
ation are two folds: First, the overall shape of ReGdi{(d) 
predicted via RG is shifted towards the smaller region 
compared to the NRG results. This leads to a smaller 
value for ojq > Q (compared to that via NRG) where 
RsGdi{\uj\ < Wo) changes its sign from positive (negative) 
to negative (positive) for < cj < cjq (— i^o < a; < 0). 
This makes the magnitudes of these additional dips and 
peaks smaller as J{uj) diverges even further (see Eq. [T8|) . 
As J is further increased, the deviations between RG 
and NRG become more transparent. This is expected as 
the perturbation theory becomes uncontrolled once the 
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FIG. 7: The Fano lineshape for pc{i^) (in unit of po) on a log- 
arithmic scale of uj/T* for different antiferromagnetic RKKY 
couplings J by NRG. The other parameters are the same as in 
Fig. 1. The dot-dash lines are power-law fits to the universal 
scaling function of lo/T* via Eq. [20l 



system moves away from the weak coupling regime. Nev- 
ertheless, the perturbative RG approach can still capture 
the qualitative features of ReGdi{i^) for T* < \uj\ < Tk 
(see Fig. [2] and Fig. [S]). Note that the perturbative RG 
approach via Eq. |6] (the Dyson's equation) can fit the 
NRG result for ReGdi{uj) better than that via Eq. [5] for 
w > Tfc, as expected. Similar to the KT scaling behavior 
for ImGdi{u!), the NRG results for ReGdi{'-^) also show 
a scaling behavior for u <T*: ReGdi{oj) « g'^g' {u /T*) 
(see Fig. |4]). Here, the scaling function g'{oj/T*) is again 
completely universal in the case of particle-hole symme- 
try. Based on the Fermi liquid theory, the very low en- 
ergy (w -C T*) crossover function for ReGdi{'-^) is linear 
in oj/T* (see, for example Eq. [9]): 



TReGdiH « -«2(^). 



(19) 



where we find 02 ~ 1.5 from the fit to the NRG result 
(see Fig. E]). 

Finally, we discuss the behavior for the Fano lineshape 
for Pc{^)- As indicated in Eq. [Tni the Fano lineshape for 
6pc{ijj) is effectively a linear combination of the asymmet- 
ric real part and symmetric imaginary part of the Gdi (w) . 
The parameter qc in Eq. [T51 depends on the conduction 
electron reservoir. Following Ref. [3| and Ref. 0, qc can 
be reasonably treated as a constant. We take a realistic 
value Qc ~ 1.4 here, corresponding to the Co/Au system 
studied in Ref. [31 and Ref. 0. We find pc{ijj) IS asym- 
metric with respect to w = with a larger magnitude 
for w > than that for w < 0. As shown in Fig. [6] and 
Fig- 13 Pc{<^) shows a dips at w « —Tk and lo T* 
as well as peaks at w ~ T/j and w k, ~T*. The peak 
(dip) at w ~ ±Tfc correspond to the first stage Kondo 



effect; while the dip (peak) at w w ±T* correspond to 
the second stage Kondo effect via RKKY coupling. We 
find a reasonably good agreement between the NRG re- 
sults and the fit via the perturbative RG approach for 
Tk ~ uj <T*. (The fit via Eq. [His somewhat better 
than that via Eq. [5] as the former gives a better fit to the 
NRG result for ReGdi{uj) The above dip-peak structure 
in the Fano lineshpe for Pc(w) in the presence of RKKY 
coupling can be detected in the STM measurement of 
the conduction electron leads as the signature of the two- 
stage Kondo effect in side-coupled double quantum dot. 
Note that the w/T* scaling in the NRG results for Pc{oj) is 
observed (see Fig. [SJ and Fig. which comes naturally 
from the scaling behaviors for both real and imaginary 
parts of Gdiii^) (see Eq. [3] and Fig. [7]). In the low en- 
ergy limit uj <^ T* where the system approaches to the 
Fermi-liquid of local spin singlet, we have the following 
approximated power- law scaling behavior for Spdi^): 



5pc{^) 
Po 



-[{l-ql)a^{ — f + 2a2q 



(20) 



We would like to make one side remark here. For 
a; ^ Tfc, the single Lorentzian approximation Eq. [S] 
can very well describe Gdi(w); however, for w < T^, we 
expect a finite contribution to 0^1(0;) from interference 
between the Kondo resonance and the broadened impu- 
rity level at dot 1. We find indeed a better agreement 
between the analytic fits and the NRG results for Gd\ (w) 
and Pciyi) via perturbative RG approach based on the 
Dyson's equation Eq. [6] than those from the single 
Lorentzian fit via Eq. [H 



V. Conclusions. 

We have studied the Fano resonance in a side-coupled 
double-quantum-dot system in the Kondo regime in the 
presence of particle-hole symmetry. In the range where 
the energy of the dot 1 is of the order of the broaden- 
ing of its energy level, quantum interference between the 
Kondo effect and the broadened energy level of the dot 1 
gives rise to modification of the Green's function on dot 
1. We apply the perturbative and numerical renormal- 
ization group approaches to describe the Fano lineshape 
in TDOS of the conduction electrons, which depend on 
both the real and imaginary parts of the Green's func- 
tion Gd\{^^ of the dot 1. At J = 0, ImGdi{oj) shows 
the Kondo peak for lo < Tk\ while ReGdi{uj) exhibits a 
peak (dip) for a; w (w w —Tk). As a result of the 
Kondo effect, the Fano lineshape in TDOS of the con- 
duction electron leads shows a peak (dip) around uj ~Tk 
(cj sa —Tk). At a finite antiferromagnetic spin exchange 
coupling between the two dots, the two-stage Kondo ef- 
fect leads to the suppression of the density of states on 
dot 1 as well as an additional dip (peak) structure in the 
real part of Gdi{uj) at w w zLT* from the NRG results. 
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This leads to an additional dip (peak) around uj ~ T* 
(w « — T*) in the conduction electron LDOS. The split- 
ing between dip and peak in LDOS at w « ±T* becomes 
more pronounced as the RKKY coupling J is increased. 
At finite values of J and for oj < T*, the NRG results for 
ReGdi{<^), ImGdiiy) and /3c(<^) all show distinct univer- 
sal scaling behaviors in lu/T*. Analytically, we find the 
perturbative RG approach can qualitatively capture the 
above behaviors for T* w < Tk- In particular, com- 
pare to the simple Lorentzian approximation for Gdi{i-^), 
we find a better fit to the NRG results for the Fano line- 
shape for /Oc(w) for T* <C a; < by taking into account 
the interference between the Kondo resonance and the 
broadened impurity level on dot 1 within the Dyson's 
equation approach. To make contact of our results in the 
experiments, the asymmetrical double dip/peak struc- 
ture and the scaling behaviors in the Fano lineshape pre- 
dicted here in the spectral properties of the TDOS of the 
conduction electron leads can be dectcted by the trans- 
port through the STM tips^i] as an indication and direct 
consequence of the two-stage Kondo effect in our side- 
coupled double-quantum-dot system. Finally, we would 
like to make a remark on the Fano lineshape in TDOS of 
the leads in our system without particle-hole symmetry. 
In this case, we expect a smooth crossover (instead of the 
KT type transition) between the Kondo and local singlet 
phases due to the potential scattering terms generated in 
the presence of particle-hole asymmetry. Nevertheless, 
further investigations via NRG are needed to clarify this 
issue. 
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